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^ : ABSTRACT 

Oh! 

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) 
observed 10 transits of XO-2b over a period of three years. We present measure- 
ments which confirm previous estimates of system parameters like the normalized 
semi-major axis (a/R*), stellar density (p*), impact parameter (6) and orbital in- 
clination (i rb)- Our errors on system parameters like a/R* and p* have improved 
by rsj 40% compared to previous best ground-based measurements. Our study 
of the transit times show no evidence for transit timing variations and we are 
able to rule out co-planar companions with masses > 0.20M e in low order mean 
motion resonance with XO-2b. We also explored the stability of the XO-2 system 
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given various orbital configurations of a hypothetical planet near the 2:1 mean 
motion resonance. We find that a wide range of orbits (including Earth-mass 



perturbers) are both dynamically stable and produce observable TTVs. We find 



that up to 51% of our stable simulations show TTVs that are smaller than the 



*0 ■ typical transit timing errors (^20 sec) measured for XO-2b, and hence remain 

O 
CO 



undetectable. 

O 



Subject headings: eclipses, stars: planetary systems, planets and satellites: fun- 
damental parameters, individual:XO-2b, methods: N-body simulations 



1. Introduction 

Observational efforts on two fronts (ground-based and space-based) have revealed a 
^reat diversity in the properties of planets (and planet candidates), and posed several new 
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questions about the planet formation process. The search for planets and the careful mea- 
surement of planetary properties are the observational foundations upon which the physics 
of planet formation may be understood. The exoplanet community has co-opted the term 
'architecture' to embody several properties of planetary systems such as multiplicity and 
orbital parameters like eccentricities (e), inclinations (i rb)> semi-major axes (a) and orbital 
periods. The distributions of individual planetary properties like masses and radii (M p and 
R P) respectively) may also be included in this term. The primary goal of the Apache Point 
Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE) is to catalog transit 
lightcurves at high precision, in order to (1) measure system parameters, and (2) look for 
transit timing variations that may indicate the presence of additional planets. 

The transit technique applies to those systems where the orbital inclination of an exo- 
planet i^_close_to_90^_J^.ed^e-on) with respect to the observer's sky-plane (see discovery 



paper ICharbonneau et al.ll2000l ). In this case, the obs erver sees a u-shaped dip in the starlight 
caused when the planet eclipses the star ( jWinnll201ll ). The objective of several ground-based 
and space-based efforts focused on transit observations is to catalog and improve measure- 
ments of system parameters, which in-turn gives us an improved picture of the architecture 
of the planetary system. This process is key toward developing theories of planet formation 
that can adequately explain the origin and evolution of all planetary systems (including our 
own). 

The APOSTLE target XO-2b is a Hot- Jupiter on a 2.6 day orbit around an early type 
K dwarf (V = 11.2, LBurke et al.l 120071). The p lanet is known to have a mass of 0.555 M Jup 
and a radius of 0.992 i?j up ( jSouthworthl l2010l ) . The system is not known to have any other 
planet-mass obje cts. Transit timing measurements from the ground are con sistent with a 
linear ephemeris ( IFernandez et al.ll2009l : ISing et al.ll201ll : ICrouzet et al.ll2012l ). although the 
cited studies note that there are statistically significant deviations in the meas urement of the 



orbita l periods. A search for additional eclipses using the EPOXI mission by iBallard et al. 



( j201ll ) had inconclusive results since the eclipse was not fully sampled. Spitzer observations 
of the secondary eclipse of XO-2b show IRAC 3.6, 4.5 and 5.8 urn flu xes which are consistent 
with the presence of a temperature inversion (IMachalek et al.l 120091 ) . There have also been 
detections of optical absorbers in the planetary atmosphere such as potassium from narrow 
band optical transmission spectrophotometry ( Sing et al.ll201ll ). Early theoretical studies 
indicated that stellar insolation levels directly influenced the presence or absence of a thermal 
inversion layer (depending on the surv ival of atmospheric absorbers iHubeny et al.l 120031 : 



Burrows et al. 



2007 1: iFortnev et al.l 120081 ). The planetary atmosphere classification system 



developed by lFortney et al.l ( 120081 ) places XO-2b in an transition zone between planets with 
(pM) and without (pL) thermal inversions. XO-2b is o ne of a handful of H ot-Jupiters 
in this region. However, observational evidence suggests ( IMachalek et al.l 120081 ) that more 



complicat ed models need to be considered. Irradiation from stellar activity may need to be 
included ( IKnutson et al.l 120101 ) and atmospheric chemist ry may also need to be considered 
to provide a more complete picture (IMadhusudhanl 120121 ) . The host star of the XO-2 system 
is known to be a high metallici ty (TFe/H] = 0.45 ± 0.02), high proper motion star (/xxot = 
157 mas/yr) in a visual binary ( IBurke et al.ll2007l ). Spectral activi ty indices show that XO-2 
is fairly inactive compared to other stars of similar spectral type ( IKnutson et al.ll2010l ). 



A primary goal of APOSTLE and other campaigns that monitor tr ansiting exoplan 
ets is to search for TTVs that reveal the presence of unseen companions ( jAgol et al.l 12005 
Holman fc Murrayl 120051 ) . In principle, Earth-mass planets in or near mean motion reso- 



nances could perturb the orbit of the transiting planet enough to produce a sinusoidal oscil- 
lation in the mid-points of transits. However, the full range of stable orbits that can produce 
a detectable TTV s ignal has never been explored. The detection of TTVs by Kepler (e.g. 
Holman et al.ll2010l ) has demonstrated that stable systems are capab le of producing TTVs, 
and other studies have explored a limited range of architectures (IHaghighipour fc Kirste 
20111 ). but the systematic exploration of parameter space of an Earth-mass planet in orbit 
near a hot Jupiter has not been undertaken. Here we examine 3.6 million possible masses 
and orbits of an approximately Earth-mass planet orbiting in or near the 2:1 outer mean 
motion resonance through N-body simulations. We find that stable orbits that also produce 
detectable TTVs do exist in the XO-2 system, and hence can exist in similar Hot- Jupiter 
systems. 

In this paper we report observations of 10 transits of XO-2b, taken as part of APOSTLE. 
In §[2] we outline our observations. In §[3] we briefly outline (i) the data reduction, photometry 
(§ 13.11) . (ii) the transit model (§ 13.21) and lightcurve fitti ng (§ 13.31) : both processes have also 
been described in previous work ( jKundurthy et al.ll2013l ). In §[4] we present our estimates of 
the system parameters for XO-2b and in the subsections § 14.11 and § 14.21 we present results 
from our study of transit depth variations (TDVs) and transit timing variations (TTVs). In 
§ [5] we present results form N-body simulations used to study the stability of hypothetical 
planetary configurations at the 2:1 mean-motions resonance. Finally, in §[6] we summarize 
our findings. 



Observations 



XO-2b was observed by members of the APOSTLE team on 10 occasions over a timespan 
of 3 years from early 2008 until the spring of 2011. All observations of XO-2 were carried out 



using Agile, a high-speed frame-transfer CCD ( IMukadam et al.ll201ll ). on the ARQj 3.5m 
telescope at Apache Point, New Mexico. The summary of observations is given in Table 
[TJ XO-2 was an early APOSTLE target and was observed using a variety of instrumental 
settings, as the team had not converged on an optim al observing; strategy prior to 2010. Early 
observations were made in the /-band (A = 805nm. ICousinslll976l : lBesselllll990l ) with several 
data sets taken with short read-out (exposure times, Column (5) in Tabled]). The short read- 
out mode allowed for fine sampling of the lightcurve, but due to the lower signal-to-noise 
and the unsuitability for characterizing systematics, this observing mode was abandoned. 
In addition, the /band images also contained a strong contribution from a fringe pattern 
due interference from the backscattering of atmospheric lines within the CCD's pixels. The 
removal of this fringe pattern is discussed in § 13.11 Early in 2010, the observing strategy 
changed to longer readouts (typically 45-75 seconds) to reduce the level of correlated noise 
(red-no ise) , and were made using the r'-band which is similar to the SDSSg t filter (A = 
626nm, iFukugita et al.lll996l ). to reduce the influence of fringes. In the long read-out mode 
the telescope was defocused to spread the stellar point-spread function (PSF) across multiple 
pixels, which minimized the systematics caused by pixel-to-pixel wandering of the PSF over 
the imperfect flatfield. The longer exposures also allowed for a greater count rate which 
maximized the signal-to-noise per image. The count rate was kept below Agile 's non-linearity 
limit of ^52k ADU and well below its saturation level of 61k ADU by small adjustments to 
the telescope's secondary focus during observations. 

XO-2 (TYC 3413-5-1) and its visual binary companion TYC 3413-210-1 (separated by 
~ 30") were the only bright stars that fit in Agile^s field of view. Both stars are of identical 
spectral type (K0V) and nearly identical brightnesses in the filters used by APOSTLE , with 
their Johnson R and I magnitudes at 10.8 and 10.5 respectively ( IMonet et al.l 120031 ). Our 
uncalibrated differential photometry also showed good agreement in their brightnesses, with 
the out-of-eclipse, un-normalized flux ratios being different by only 4% and 3% in the I- 
band and /-band, respectively (Column 'Flux Norm.' in Table [1]). The observations were 
made over a variety of observing conditions (Column 'Obs. Conditions' in Table [1]). The 
observing conditions are classified as 'Clear' or 'Poor Weather' with the former implying 
good data with few or no interruptions in data collection, and the latter indicating that the 
we experienced cloud cover or poor seeing conditions resulting in lower quality data. The 
tabulated transits are those for which we were able to capture the whole transit or at least a 
partial transit. Partial transits are those where portions of the in-eclipse lightcurve were lost 
due to bad weather or instrumental failure. Several data points were lost for the nights of 
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UTD 2008-09-22 (#4), 2009-02-06 (#5), 2009-03-12 (#6) and 2010-12-27 (#8). However we 
do include these nights since we did manage to obtain reasonable portions of the in-eclipse 
and out-of-eclipse data, which make it possible to determine transit properties (albeit at a 
loss of accuracy and precision) . 



3. Data Analysis 

This section outlines various stages in the analysis of lightcurves, starting with (i) the 
image reduction a nd photometry, and (ii) the transit model and Markov Chain Monte Carlo 



analyzer (see also. iKundurthy et al.ll2013l ). 



3.1. Reduction & Photometry 

APOSTLE data were reduced using a pipeline developed specifically for Agile images. 
The pipeline (written in IDlJj) performs pixel-by-pixel error propagation, and image pro- 
cessing specific to the Agile CCD. In addition to the standard photometric reduction steps 
like dark subtraction and flat-fielding, the pipeline performs non-linearity and fringe cor- 
recti ons specific to Ag i le. Th e details on Agile' s non-linearity correction are described in § 



3 of IKundurthy et al.l ( 120131 ). Some of the initial observations by APOSTLE were carried 
out using short exposures (see Column (5) 'Exp.') in Table [H We binned these data by 
averaging the flux ratios in 45 sec bins. 

For several of the initial XO-2 data taken using Agile' s /-band, a fringe pattern had 
to be subtracted to create science images. Photons from strong atmospheric lines (in the 
/-band bandpass) backscatter within the CCD pixels, and owing to the variable thickness 
of the pixel array, interference between these photons creates the fringe pattern. Since the 
fringe pattern is convolved with the illumination pattern of the CCD, the fringe correction 
has to be applied after dark-subtraction and flat-fielding. During the initial characterization 
of the Agile CCD, an empirical fringe frame was produced by median combining dithered 
frames on the dark sky, where the flux contributions from stars were removed by outlier 
rejection. The resulting combined frame served as a "Model" fringe pattern (F), normalized 
to have a median of zero, and an amplitude of one, such that it could be scaled to match the 
fringe patterns on science frames. The science frame affected by the fringe pattern (T ) is 
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assumed to be a linear combination of the fringe-less science image (T) and a fringe pattern: 

T f = T + a F (1) 

where F is the model fringe frame and a is a scaling factor describing the amplitude of the 
fringing on a given frame. The fringe amplitude is estimated by minimizing \ 2 an d fitting for 
a in the above model. The corrected images (T = T — clqF), were found to be sufficiently 
corrected of fringes after visual examination. We found the fringe amplitudes (ag) on XO-2 
science frames to always be smaller than the standard deviation (i.e. scatter) of the global 
sky background on each frame; a ranged between 6-13% of the scatter in the sky for XO-2 
data. 

We extracted photometry from an optimal circular aperture centered around the tar- 
get and comparison stars. In addition we extracted the counts on image products like 
the master-dark and master-flat, to serve as nuisance parameters for detrending, using the 
same aperture and centroids from photometry. The centroids were derived using SExtractor 



(IBertin fc Arnoutsl Il996l ) , which allowed for the use of customized donut-shaped convolu- 
tion kernels for defocused PSFs. Science frames where pixels inside a photometric aperture 
exceeded Agile' s saturation limit of 61k were rejected. Images at the other extreme, where 
the stars were obscured by clouds, and resulted in low signal to noise measurements were 
also rejected (i.e. where individual photometric errors were > 5000 ppm). The fraction of 
rejected frames per night is listed in Column (9) '%Rej.' in Table [TJ The optimal aperture 
was selected after extracting photometry on a list of circular apertures with radii between 
5-50 pixels at an interval of 1 pixel. The optimal aperture was selected where the scatter 
in the residuals of the detrended lightcurve minus a trial transit model (based on values 
from the literature) was minimized. The correction function F cor (or detrending function) is 
modeled as a linear sum of nuisance parameters as described by the following equation: 

" cor,i / J ^k^k^ii \^J 

fc=l 

where X^i are the nuisance parameters, c^ are the corresponding coefficients. The index 
k counts over the number of nuisance parameters N nus . The detrending coefficients are 
chosen by minimizing the x 2 between the observed data (O), a model function (M) and the 
correction function, 



x 2 = ^ (Oj ~ Ms- F covJ ) (3) 

j i 

here j is the index which counts over the total number of data points (iV a ii). The observed, 
model and correction function terms are all in normalized flux ratio units. The list of param- 
eters used for detrending each lightcurve in the XO-2 dataset are presented in Table El with 



variable definitions given in the footnotes. A suitable set of detrending parameters were se- 
lected for a given night by running several manual trials of linear least squares minimization 
on single transit lightcurves and its corresponding model lightcurve and correction function 
(Eq. [2j). The model parameters were fixed at reasonable values and only the coefficients 
(cfc) were fit. Those nuisance parameters which returned large uncertainties on the coeffi- 
cients were excluded since this indicated a poor match to the noise trend in a lightcurve. 
This method of selecting nuisance parameters is admittedly ad-hoc; automated parameter 
selection techniques could be applied in the future. The detrending is performed using the 
final set of nuisance parameters in conjunction with fitting the transit parameters in order 
to ensure that the correction process is not biased by the trial model used in the selection 
of nuisance parameters. 



Table 1: APOSTLE Observing Summary for X02 



T# 


UTD 


Obs. Cond. 


Filter 


Exp. 


Bin 


Phot. Ap. 


RMS (ppm) 


%Rej. 


Flux Norm. 


Error Scaling 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


2008-01-09 


Clear 




0.5 


45 


13 


557 


3% 


0.9561 


1.1740 


2 


2008-02-12 


Clear 




0.5 


45 


24 


510 


<1% 


0.9561 


1.3129 


3 


2008-03-04 


Clear 




0.5 


45 


15 


411 


<1% 


0.9632 


1.2510 


4 


2008-11-23 


Poor Weather 




10,25,45 


45 


27 


939 


1% 


0.9756 


2.7482 


5 


2009-02-07 


Poor Weather 




45 


- 


31 


1085 


11% 


0.9721 


2.9868 do 


6 


2009-03-13 


Poor Weather 




45 


- 


46 


405 


11% 


0.9766 


1.2356 ' 


7 


2010-10-25 


Clear 


r' 


45 


- 


35 


553 


2% 


0.9676 


2.2760 


8 


2010-12-27 


Poor Weather 


r' 


45 


- 


43 


775 


5% 


0.9643 


3.0092 


9 


2011-01-30 


Clear 


r' 


45 


- 


48 


354 


1% 


0.9646 


1.4571 


10 


2011-03-05 


Clear 


r' 


45 


- 


43 


693 


<1% 


0.9660 


2.7990 



(I) Transit Number, (2) Universal Time Date, (3) Observing Conditions, (4) Observing Filter, (5) Exposure Time (seconds) 
(6) Bin size in seconds, (7) Optimal Aperture Radius (pixels), (8) Scatter in the residuals 

(9) % frames rejected due to saturation or other effects, (10) Flux normalization between the target and comparison star 

(II) The factor by which the photometric errors were scaled 



Table 2: Nuisance Parameters used for XO-2 Nights 

~T# Nuisance Parameters used (FLDC, OLDC and MDFLDC) 



Tl 
T2 
T3 
T4 
T5 
T6 
T7 
T8 
T9 



airmass, mskyi, msky 2 , gsky, xi, yi, x 2 , y 2 , sDi, sD 2 , sFi, sF 2 , a 
airmass, mskyi, msky 2 , gsky, xi, yi, x 2 , y 2 , sDi, sD 2 , sF x , sF 2 , a 
airmass, mskyi, msky 2 , gsky, x u y u x 2 , y 2 , sDi, sD 2 , sF x , sF 2 , a 
airmass, mskyi, msky 2 , gsky, x b y b x 2 , y 2 , sDi, sD 2 , sF x , sF 2 , a 
airmass, xi, y 1? sDi, sFi, ao 
airmass, xi, y 1? sDi, sFi, a 
airmass, gsky, xi, y 1? sDi, sFi 
airmass, gsky, xi, y 1? sDi, sFi 
airmass, gsky, xi, yi, sDi, sFi 



T10 airmass, gsky, xi, y 1? sDi, sFi 



airmass = atmospheric column 

( x i 5 yi)? ( x 2 5 J2) — Centriods of Target (1) and Comparison (2) 

mskyi, msky 2 = Median sky around Target (1) and Comparison (2) 

gsky = Global median sky 

sDi, sD 2 = Sum of counts in aperture on Master-Dark 

sFi, sF 2 = Sum of counts in aperture on Master-Flat 

ao = Fringe scaling (see Eq. []]) 
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Table 3: APOSTLE Lightcurve data* for XO-2 



T# 


T-TO 


Norm. Fl. Ratio 


Err. Norm. Fl. Ratio 


Model Data 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


-0.1341616 


1.000004100 


0.001207284 


1.000000000 


1 


-0.1315574 


0.999797978 


0.001102731 


1.000000000 


1 


-0.1305157 


0.997707769 


0.001047412 


1.000000000 


1 


-0.1299949 


0.999906661 


0.001190794 


1.000000000 


1 


-0.1289532 


0.998693635 


0.001118320 


1.000000000 


1 


-0.1279116 


0.998939727 


0.001098241 


1.000000000 


1 


-0.1273907 


0.999021973 


0.000987372 


1.000000000 


1 


-0.1268699 


1.001381660 


0.001209313 


1.000000000 


1 


-0.1263491 


0.998899046 


0.001076096 


1.000000000 


1 


-0.1258282 


1.001212353 


0.001160014 


1.000000000 



*The data are presented in their entirety as an online-only table. 
(1) Transit Number, (2) Time Stamps - Mid Transit Times (BJD) 
(3) Normalized Flux Ratio, (4) Error on Normalized Flux Ratio 
(5) Model Data 
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The 10 transits of XO-2b are shown in Figure [T] in normalized flux ratios (with offsets 
for clarity). The plotted data result from the data reduction and model fitting processes 
described in sections § 13. 11 § 13.21 and § 13.31 



3.2. MultiTransitQuick 



We developed a transit model called MultiTransit Quick (MTQ) in 



is based on the analytic lightcurve models presented in iMandel fc Agol 



PYTHON implementation of some of its functions (from EXOFAST bv lEastman et al.ll2012l). 



'YTHON, which 



(120021). and th e 



The d escription of MTQ used for this study is described in some detail in iKundurthy et al. 



( 120 13l ). MultiTransitQuick can be used to (1) fit for transit depths for data taken with 
multiple filters (Multi-Filter model), or (2) fit each transit depth individually (Multi-Depth). 

The set of parameters used for Multi-Filter version of MTQ is #Muiti-Fiiter = {^r, tc ) 
Dj...n f , ^ij...Nn V2j...n f , Ti...Ar T }, where t F is transit duration, to is the limb-crossing duration 
and Ti are the mid-transit times. The filter dependent parameters include the transi t dept h 
Dj and the limb- darkening parameters vij and v 2 j described in IKundurthy et al.l (120131 ). 
The subscripts i...N T and j...N F are used to denote multiple transits (N T ) and multiple 
filters (N F ) respectively. The XO-2 data were gathered using the /-band and /-band, 
where N F was 2 and the number of transits N T was 10. 

The parameter set is only slightly different for the Multi-Depth version, with ^Multi-Depth 
= {^t, tG-> D i _ NT) Vij...n F i V2j...n f , T1...at t }; the difference being the transit depth is now fit 
for each transit separately instead of each filter. The filters are still tracked to ensure the 
use of the correct limb- darkening parameters for each lightcurve. A single lightcurve is set 
as the "reference" lightcurve and used to internally compute several orbital parameters. For 
the XO-2 data set we used # 3 since it had few gaps, and had the best photometric precision 
among the /-band lightcurves. We did not use an /-band lightcurve as a reference since this 
filter is more affected by limb-darkening when compared to the /-band and hence estimates 
of the planet-to-star size ratio (R p / R*) (used in the computation of orbital parameters) may 
be susceptible to degeneracies. Using the Multi-Depth model aids in understanding transit 
depth variations over multiple epochs. 



3.3. Transit MCMC 



We developed a Markov Chain Monte Carlo (MCMC) analyzer called Transit MCMC 
(TMCMC), based on the Metropolis-Hastings (M-H) algorithm (IGelman et al.ll2003l : lTegmark et al 
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2004J : lFordl[2005h . with an adaptive step-size modifier (jCollier Cameron et al.1l2007h . MCMC 
routines quantify the uncertainty distributions of model parameters (given the data) using 
Bayes' theorem, by sampling parameter space such that samples from high probability re- 
gions (low x 2 ) are selected at a greater rate than those from low probability regions. The 
final ensemble of sampled points the MCMC routine typically represent the uncertainty dis- 
tributions of the model parameters. One must note that adaptive MCMCs are generally not 
considered to be trul y Markovian in nature and their results are valid only if adaptation 
diminishes with time ( [Roberts fc Rosenthal! 120091 ): a property that our chains do display. 



For APOSTLE data sets, we explored system parameters using three different kinds 
of chains. Two of these were based on the Multi-Filter parameter set #Muiti-Fiiter described 
in § 13.21 and third used the Multi-Depth parameter set (^Multi-Depth)- The two Multi-Filter 
chains used Fixed Limb- Darkening Coefficients (LDC) and Open Limb- Darkening Coeffi- 
cients (OLDC). For the Fixed LDC chains (FLDC ), the coefficients were simply fixed to 
values tabulated for the appropriate observing filter ( IClaret fc Bloemenll201lh . For the Open 
LDC chains, the limb- darkening parameters v± and i?2, for both /-band and riband data, 
are allowed to float. The ability to constrain stellar limb-darkening: requires high pr ecision 
data, such as those collected using the HST ( IBrown et al.ll200ll : iKnutson et al.l 120071 ). Pre- 
vious attempts to fit for limb-darkening on APOSTL E data have resulted in Markov chains 



that failed to converge (jKundurthy et al. 



2011 



20131 ). The third type of Markov chain was 



run on the Multi-Depth parameter set ^Multi-Depth described in § 13.21 APOSTLE lightcurves 
were gathered over a long time-baseline, and statistically significant depth variations seen 
in the data may help shed light on the various phenomena responsible for depth variations 
(see § 13.21) . or point to limitations in the data and model. 

Several of the prel iminary steps for executing a chain using TMCMC are described in 
( JKundurthy et al.l 120131 ). These steps include (1) setting bounds and (2) running short 
single-parameter exploratory chains to determine a set of suitable starting jump-sizes for 
the long Markov Chains. We ran long chains of 2 x 10 6 steps from two different starting 
locations for each model scenario: Fixed LDC, Open LDC and Multi-Depth/Fixed LDC. 
After completion we (1) cropped the initial stages of these chains to remove the burn-in 
phase, where the chain is far from the best-fit region, and (2) we exclude the stage where 
the chai n is far from the op timal acceptance rate of 23± 5% (as noted for multi-parameter 
chains, iGelman et al.l 120031 ). We run three types of post-processing on the chains after 
cropping: (a) We compute the ranked and unranked correlations in the chains of every fit 
parameter with respect to the others. These statistics provide an estimate of the level of 
degeneracy between parameters in a given model. The next post-processing steps are two 
commonly used diagnostics to check for chain convergence, na mely (b) computing the auto- 
correlation lengths and (c) the Gelman- Rubin R-static values ( IGelman fc Rubinlll992l ). The 
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auto-correlation lengths determine the scale over which a chain has local trends. From the 
auto-correlation length one can compute the effective length as the total chain length divided 
by the auto-correlation length, which represents the statistical significance with which the 
uncertainty distribution was sampled. A large effective length (> 1000) represents a well- 
sampled distribution. The R-statistic represents the level of coverage the chain has over 
the parameter space. When parameter space has been properly sampled the R-statistic 
computed using chains from different starting locations will be close to 1. We deem those 
chains as coverged that have an effective length > 1000 and an R-statistic within 10% of 1. 



3.3.1. TAP 



The Transit Analysi s Package (TAP: lGazak et al.ll201ll ) implements the red-noise model 



of ICarter fc Wind (120091 ) . who find that models that do not fit for red-noise are subject to 
inaccuracies in transit parameters on the order of 2-3a and tend to have underestimated 
errors by up to 30%. For transit timing studies, poor estimates such as these are cause for 
concern, since smaller errors and large deviations from the expected time can easily lead 
to false claims of TTVs. Since TMCMC does not include red-noise analysis we run fits on 
APOSTLE lightcurves using TAP as a check to the results derived from TMCMC. 

The typical TAP parameter set is: T ap = {a/R*, i or b, (Rp/R*)i...N F , T immmNT , cr( w hite,;...iV T ), 
^(red,i..JV T )}j where a/R*, i or b and (R p / R*) are the commonly fit transit parameters denoting 
the semi-major axis (in stellar radius units), the orbital inclination and the planet-to-star 
radius ratio, respectively. The noise analysis parameters (J( w hite,i...N T ) an d &(red,i...N T ) are the 
white-noise and red-noise levels for N T transits respectively. The TAP package does not fit 
for the period using the transit times, and often yields poor estimates of the period, so we 
fixed the period to the value derived from TMCMC. The limb- darkening was fixed to values 
from the literature. The orbital eccentricity (e) and argument of periastron were kept fixed 
at for XO-2b. 



4. System Parameters 

This section describes results from our execution of the 2 chains for the parameter 
set ^Multi-Filter, and 1 for the ^Multi-Depth described in § 13.31 Post processing statistics and 
other data for these chains are listed in Tabled! The columns Wf ree ', 'Chain Length', 'Corr. 
Length' and 'Eff. Length' list the number of free parameters, the length of the cropped chain, 
the correlation and effective lengths, respectively. All chains were run for approximately 2 
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Million steps, but about 100,000 of the initial steps were removed to account for "burn- 
in" and selection rate stabilization. The XO-2b chain with Open LDC has a low effective 
length indicating poor Markov chain statistics. The Fixed LDC chains (both 0Muiti-Fiiter an d 
^Multi-Depth) model satisfy the condition of a well sampled posterior distribution (effective 
length is > 1000). The final two columns list the goodness of fit (i.e. lowest x 2 m the 
MCMC ensemble) and Degrees-of-freedom (DOF) from the respective chain. Parameters 
from all chains had Gelman-Rubin R-statistics close to 1 indicating that the parameter 
space was covered evenly (though the OLDC chain was not sampled finely enough, based on 
the auto-correlation data). 



Table 4: TMCMC Chains for XO-2 



Chain 


Model Vector 


Nfree 


Chain Length 


Corr. Length 


Eff Length 


x 2 


DOF 


FLDC 

OLDC 

MDFLDC 


^Multi-Filter 
^Multi-Filter 
^Multi-Depth 


14 
18 
22 


1,900,001 
1,900,001 
1,900,001 


100 
4,949 

588 


19,000 

383 

3,231 


3272.94 
3268.48 
3197.53 


3312 
3308 
3304 



Table 5: XO-2 Parameters for 0Muiti-Fiiter 



Parameter 


FLDC 


OLDC 


Unit 


MTQ Parameters 


to 


0.0107±0.0002 


0.0108±0.0002 


days 


t T 


0.1008±0.0001 


0.1004±0.0002 


days 


D (D 


0.0126±0.0001 


0.0126±0.0001 


- 


2? (r .) 


0.0130±0.0001 


0.0131±0.0001 


- 


Mi) 


(0.5944) 


0.5440±0.0313 


- 


vi (r') 


(0.6980) 


6411+ 00332 


- 


Mi) 


(0.1452) 


0.2806±0.0904 


- 


v 2 (r') 


(0.3524) 


0.478618:8^ 


- 


Derived Parameters 


(Rp/R*)^ 


0.1030±0.0003 


0.1033±0.0004 


- 


(Rp/R*)^ 


0.1024±0.0003 


0.1029±0.0004 


- 


b 


1 7 +0 - 04 


0.21±0.03 


- 


a/R* 


8.14±0.06 


8.11±0.07 


- 


i or b 


88.79±0.15 


88.531™ 


°(deg) 


v/R* 


19.55±0.15 


19.48±0.17 


days -1 


P* 


1.49±0.03 


1.48±0.04 


g/cc 


P (2.6159 days +) 


-3467±22 


-3466±22 


milli-sec 
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Table 6: XO-2 Parameters for ^Multi-Depth 



Transit Depths 



Value 



Units 



Rp/R* 



Value 



Units 



/-band 



(D)i 
(Dh 
(D)s 
(D)a 
(D) 5 
(D) e 



0.01338±0.00013 
0.01300±0.00017 
0.01226±0.00013 
0.01199±0.00013 
0.01209±0.00039 
0.01267±0.00015 



(Rp/R*) i 
{R P /R*)2 
(Rp/R*) 3 
(Rp/R*)^ 

(R P /R*) 5 
(R P /R*)6 



0.1061±0.0005 
0.1047±0.0007 
0.1019±0.0005 
0.1009±0.0005 
0.1013±0.0015 
0.1035±0.0006 



-band 



(D) 7 
(D)s 
(Dh 
(D)io 



0.01326±0.00009 
0.01320±0.00015 
0.01307±0.00009 
0.01300±0.00013 



(Rp/R*)7 
(R P /R*)8 

(Rp/R*)g 

(R P /R*)io 



0.1032±0.0004 
0.1030±0.0006 
0.1026±0.0003 
0.1024±0.0005 



Other MTQ Parameters 


Parameter 


Value 


Units 


Parameter 


Value 


Units 


t G 

vi(l) 

V2(l) 


0.0108±0.0001 

(0.5944) 
(0.1452) 


days 


t T 

^i(r') 


0.1008±0.0001 
(0.6980) 
(0.3524) 


days 


Derived Parameters 


b 

P* 


0.24±0.02 

88.29±0.15 

1.43±0.02 


°(deg) 

g/cc 


a/R* 
vjR* 


o n9 +0.03 
iq27 +0.08 

iy.z/_ 10 


days 1 



The resulting best-fit parameter estimates are listed in Table [5] for the Multi-Filter 
models, and in Table [6] for the Multi-Depth models. These tables also list the derived 
system parameters. The tr ansformation betwe en th e MTQ parameters to the derived system 
parameters are described in lCarter et al.l ( 120081 ) and iKundurthy et al.l ( J201 lh . Contour plots 
showing the joint probability distributions (JPDs) for the fit and derived parameters are 
shown in Figures [2] and [3l respectively. There are no strong correlations between the fit 
parameters £j-, tG an d the transit depths, Dj and D r /, as seen in Figured 
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APOSTLE Observations of X02 




-0.150 -0.075 0.000 0.075 

BJD - Mid Transit Time (days) 



0.150 



Fig. 1. — Six /-band and four r'-band lightcurves of XO-2b. The vertical axis is in nor- 
malized flux ratio units. The horizontal axis shows time from the mid-transit time in days, 
computed by subtracting the appropriate mid-transit time for each transit from the best-fit 
values in the Fixed LDC chain. 
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XO-2 FLDC 



u 1.50e+04 
| 1.00e+04 
Z 5.00e+03 



0.1009 






0.1005 



0.0127 



0.0124 



0.0131 



0.0129 




0.0106 0.0111 

to 



0.1005 0.1009 0.0124 0.0127 



1.50e+04 
Number 



Fig. 2. — Plots of the Joint Probability Distributions (JPD) of parameters from the Fixed 
LDC chains, demonstrating that parameters chosen for #Muiti-Fiiter are generally uncorrelated. 
Table [5] gives units. 
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JPD - XO-2 System Parameters 



g 1.50e+04 
'I 1.00e+04 
£ 5.00e+03 



Fernandez et al. (2009) 
Crouzet et al. (2012) 
Sing et al. (2011) 
APOSTLE (TMCMC) 
Burke et al. (2007) 




1.50e+04 £ 
1.00e+04 1 
5.00e+03 Z 



.00 8.19 88.40 89.45 0^8 0.22 L42 U>2 
a/R* iorb b P* 



1.50e+04 
Number 



Fig. 3. — Plots of the Joint Probability Distributions (JPD) of derived system parameters 
from the Fixed LDC chains. Parameter estimates available in the literature are overplotted. 
Table [5] gives units. 



System parameters agree with previously published values in the literature, as seen 
by the overlap of the uncertainties in the JPD plot (Figure BJ). In addition APOSTLE's 
measurements give tighter constraints on several of the system parameters. The errors 
on a/R*, i or b, the impact parameter (b) and the stellar density (p+) are more precise than 
previ ous measurements by a factor of ^3 ( IBurke et al.ll2007l ; iFernandez et al.ll2009l ; ISing et al. 



201 lh . However, previous studies using the combination of the TMCMC Markov chain analyzer 
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and MTQ t ransit model have resu lted in underestimated errors since we do not include rednoise 
analysis ( jCarter fc Winnll2009l ). Hence, the errors presented in Table [6] and Figure [3] do not 
truly reflect improvements in the measurements of parameters. More conservative constraints 
were placed on a subset of these system parameters using the TAP package. Comparisons 
of some parameters and their uncertainties are presented in 



Table 7: Comparison of Estimates of System Parameters for XO-2b 



Parameter TMCMC TAP 


B07 


F09 Sll 


C12 Units 


a/R* 8.14±0.06 8.13±0.10 


8.20±0.20 


8.13±0.20 7.83±0.17 


7.99±0.07 


i orb 88.79±0.15 88.80±0.61 


88.90±0.75 


87.62±0.51 


88.01±0.33 °(deg) 


b 0.172to:o2i 0.171±0.085 


0.158±0.110 


0.160±0.110 0.324±0.070 


0.280±0.044 


p* 1.49±0.03 1.49±0.05 


1.52±0.11 


1.48±0.10 1.33±0.09 


1.41±0.04 g/cc 


TMCMC & TAP values are from independent analysis of APOSTLE lightcurves 

B07 - (Burke et al. 2007). F09 - (Fernandez et al. 2009). Sll - (Sine et al. 2011). C12 - 


fCrouzet et al. 2012) 



It is clear that using TAP on the APOSTLE dataset and accounting for rednoise provides 
more conservative estimates of the system parameters when c ompared t o TMCMC values. 
TAP errors for a/R* and p* are better than those reported by ISing et al.l (1201 ll ) by upto 
^■40%; whose obse rvations we r e mad e using the 10.4m Gran Telescopio Canarias (GTC). 
One must note the ISing et al.l ( 120111 ) data result from narrow-band photometry and hence 
would have much lower photometric precision when compared to broadband observations 
from the same telescope. The resulting TAP errors for the impact parameter (6) and the 
orbital inclination are larger than those from the GTC study by ~ 20%. Thus we report 
only some improvements in the measurements of sy stem paramete r s. Ou r system parameters 
agree with the current best estimates reported in ICrouzet et al.l ( 120121 ) , but our error bars 
are larger by factors < 2. 



4.1. Transit Depth Analysis 

Figure [4] shows transit depth vs. transit epoch for the /band (top panel) and riband 
(bottom panel). The overall variations in the /-band depth are ^ 0.05% compared to the 
0.01% uncertainty in D(j) from the joint fit to depths in the Fixed LDC chain. Depth 
variations can be caused by spots. Even though the variations we present are significant we 
refrain from claiming the detection of spot-modulation. These deviations are likely due to 
the incomplete sampling of several transits. Spots influence stellar brightness to a greater 
extent at shorter wavelengths, so the r'-band would be more conducive to showing depth 
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Transit Depth Variations XO-2 



• MDFLD values 
Weighted Mean of MDFLD values 

Depth from FLD 

— — 1-cr errors from FLD 




0.0120 



0.0135 
0.0134 
0.0133 
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0.0131 
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Fig. 4. — The transit depth D as a function of transit epoch for both /band and r'-band 
observations of XO-2b. The solid horizontal and dashed lines represent the best-fit value 
and errors respectively for D from the Fixed LDC TMCMC fit. The dotted line is the weighted 
mean of transit depth values from the Multi-Depth Fixed LDC chains. 



variations. However, the overall depth variation in the /-band lightcurves is of the order 
0.01% and is consistent with the errors on D( r /) from the fit reported by the Fixed LDC 
chain (Table Ej). The variations seen in the /-band depths are difficult to exaplin. They can 
either be due to (1) real brightness variations caused by spot-modulations, or (2) variations 
arising from the transit model's inability to accurately constrain transit depths and errors 
given the incomplete sampling of lightcurves (like Transit # 4, 5 and 6). We note that 
errors on the transit depth from our Multi-Depth LDC chains are more sensitive to the 
incompleteness of in-eclipse data rather than the photometric precision of a given lightcurve. 
For example, transit numbers #4 and #6 have 86% and 72% of the in-eclipse transit data 
sampled respectively; they also have lightcurve residuals of 939ppm and 405ppm respectively. 
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So even though # 6 has significantly better photometric precision, the depth for # 4 is 
constrained slightly better due to the fact that more in-eclipse points are available. A 
more dramatic difference can be seen with transit #5, which has only 47% of the in-eclipse 
lightcurve sampled and the worst depth constraint. Other than #5 all transits have eclipses 
sampled to better than 70%, hence the depth uncertainties are all comparable. Even though 
the Multi-Depth LDC chain satisified the convergence criteria described in § 13.31 the effective 
length of the chain is far lesser than the value derived for the Multi-Filter LDC chain (see 
Table Hj), indicating the lower significance of these MCMC results. In addition one must 
note the lack of rednoise analysis in TMCMC, which implies that the errors presented Figure [4] 
are probably underestimates. 



4.2. Transit Timing Analysis 



Several planetar y systems have been ob s erved to have Transi t Timing; Variations (T TVs) 
faolman et all I2OI0I : llissauer et all boilal : ballard et all l201ll : iNesvornv et all I2OI2I ). In 
certain configurations these variations can be on the order of minutes and can be easily seen 
in an diagram showing observed minus computed transit times (i.e. the difference between 
the measured times and times expected from purel y Keplerian orbita l per iods) . Using TTVs 
to loo k for additional planets was first proposed by lAgol et al.l ( 120051 ) and lHolman fc Murray 
( 120051 ) , who showed that unseen planetary siblings can gravitationally influence the orbits of 
a known transiting planet. The TTV signals are known to be especially strong if the unseen 
companion lies close to mean motion resonance with the transiting planet. 

Available transit times from the literature for XO-2b include those of the disc o very paper 



bvlBurke et al.l (2007). and the f ollow-up observations by Fernandez et al J j 20091). ISing et a 



(12007 



( j201ll ) and lCrouzet et al.l ( 120121 ). We excluded the 11 transit times from lBurke et al 
since the reported timing errors were on the order of ~ 3min, and are too large to provide 
meaningful constraints on the orbital ephemeris. We do include transit times from the 
other three studies. The ti me coordinate BJD (T DB) has become the standard system used 
for transit timing studies ( jEastman et al.l l2010l ) and all transit times used in this study 
were brought to this system. The timestamps of all APOSTLE data w ere converted to 
BJD (TDB) in the customized reduction pipeline (JKundurthv et al.lboill ). The APOSTLE 
pipeline's time conversions have been verified by comparis on to the commonly used time 
conversion routines made available by Eastman et al.l ( 120 lOl ) . 



The observed minus computed mid-transit times (O-C) are plotted in Figure [51 There 
are two versions of APOSTLE transit times presented, one from the Fixed LDC chain (TMCMC 
+ MTQ), and the other from the TAP fit to APOSTLE lightcurves. These transit times are 
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Table 8: APOSTLE Transit Times for X02 



Epoch TO (TMCMC) 


cro 


TO (TAP) 


<7xo 




2,400,000+ (BJD) 


(BJD) 


2,400,000+ (BJD) 


(BJD) 


26 


54474.73242 


0.00011 


54474.73252 


0.00021 


39 


54508.73875 


0.00013 


54508.73877 


0.00020 


47 


54529.66596 


0.00009 


54529.66604 


0.00016 


148 


54793.86845 


0.00024 


54793.86836 


0.00041 


177 


54869.72753 


0.00023 


54869.72735 


0.00050 


190 


54903.73397 


0.00014 


54903.73403 


0.00020 


416 


55494.91796 


0.00015 


55494.91794 


0.00021 


440 


55557.69891 


0.00020 


55557.69892 


0.00025 


453 


55591.70507 


0.00008 


55591.70511 


0.00015 


466 


55625.71093 


0.00014 


55625.71112 


0.00019 


Fit 


3 eriod (days) 


<7p 


TO (BJD) 


CT0 


TMCMC 


2.615860095 ± 0.000000209 


2454474.7327333 


± 0.0000599 


TAP 


2.615860014 ± 0.000000346 


2454474.7327964 


± 0.0001028 



tabulated in Table El A linear ephemeris was fit to all the data (including literature values) 
using the equation, 

Ti = TO + Epoch- x P (4) 

resulting in a best fit ephemeris of, 

P = 2.61585988 ± 0.00000016 days 
TO = 2454406.720516 ± 0.000046 BJD 



with a goodness of fit x 2 = 105.14 for 31 Degrees-of- freedom. The reduced chi-squared 
X 2 /v = 3.39 indicates that the linear fit is not robust, and is either due to the l arge timing 



devia tions in the data or underestimated errors. For example, several of the ISing et al. 



( j201ll ) data points lie far from the zero O-C line, with the largest deviation being ^ 106 
sec. The scatter in the O-C values as a whole is <~ 39 sec which makes the 106 sec deviation 
fall within the 3a confidence interval of the collective data set. Though the linear ephemeris 
does not precisely fit the timing data, the level of variation is not significant enough for us 
to claim unseen planets as the cause. 

In order to compare the ephemerides derived with and without rednoise analysis, we 
fit for a linear ephemeris to the APOSTLE transit times from TMCMC and TAP respectively 



150 



100 



50 



23 



Transit Timing Residuals for XO-2 
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Fig. 5. — The Observed minus Computed Transit Times for XO-2b. Values from APOSTLE's 
TMCMC fit, TAP and the literature are plotted. The horizontal axis represents the transit 
Epoch. The zero-line ephemeris is described in § 14.21 



(presented in Table [8j). The difference between the periods derived for these subsets and the 
period derived from all available transit times was < 7 milli-seconds. The results from the 
fits to the TMCMC and TAP transit times are presented at the bottom of Table [51 The reduced 
X 2 s for linear ephemeris fits to the TMCMC and TAP subsets were 3.96 and 1.26 respectively, 
confirming that TAP gives more conservative errors for the transit times thanks to the red- 
noise analysis. 

We were also able to rule out sinu soidal trends in the data by running a generalized 
lomb-scargle analysis on the O-C data ( Zechmeister fc Kiirsterl I2009J ) . We fit for a period, 
amplitude and phase offset on two sets of O-C data. The first set included the literature 
dataset and all APOSTLE measurements, and the second included only the TAP measure- 
ments of APOSTLE. We found that a sinusoid of period <~ 19 days and amplitude of <~ 32 
seconds improved the fit when the entire timing data set was used. The sinusoid fit yielded 
a x 2 =81.8, which was an improvement of Ax 2 = 23.3 compared to the linear ephemeris (see 
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above). However, the reduced x 2 remains greater than 1, making it a non-robust fit. The 
periodogram analysis on only the TAP data, yielded a period and amplitude of 8.5 days and 
27.4 seconds respectively. The A% 2 sh owed improvement compared to the linear ephemeris 



model by ^7.6. In a manner similar to lBecker et al.l ( 120131 ) we tested the significance of this 



sinusoidal fit by repeating the analysis on 10 5 randomly cycled permutations of the ampli- 
tudes of the O-C measurements, keeping the epochs fixed. We found that 76.7% of the fits 
showed improved A\ 2 at amplitudes greater than or equal to 27.4 seconds, indicating that 
the periodicity is not likely from a real TTV. 



For the case when planets are in mean-motion resonance (MMR), lAgol et al.l (120051 ) 



showed that the analytic expression, 5i max *~ -^ (m m I^ y, can roughly estimate the 

amplitude of the timing deviation (6t max ). The quantities m pert , m tra ns, P and j are the 
mass of the unseen perturber, the mass of the transiting planet, the orbital period of the 
transiting planet and the order of the resonance respectively. For the XO-2 system, we can 
rule out possible system configurations given the amplitude of the weak sinusoidal fit to the 
TTV data (^ 27.4 sec). Using the orbita l period from Table M and the mass of XO-2b, 



^trans = 0.565 M Jup ( IFemandez et al.ll2009l ). we compute the maximum mass perturber that 
could exist in the XO-2 system in the 2:1 MMR to be ~ 0.2M e , i.e. additional planets with 
M p < 0.2 M e may exist near the 2:1 MMR given these data. At higher order resonances, 
this maximum mass (for a possible perturber) is larger. A more detailed analysis of TTV 
needs to account for a variety of orbital configurations of hypothetical companions, which is 
addressed in the following section. 



5. System Dynamics 

In order to evaluate the stability of systems that produce a detectable TTV, we in- 
tegrated 3.6 x 10 6 orbital configurations designed to mimic the XO-2 system. Each trial 
consisted of the known planet XO-2b and a hypothetical terrestrial-like companion. 

The parameter space we cover is presented in Table [91 In this table A is the interval 
between values, all of which were varied uniformly in the range between "Min" and "Max," 
producing N bins. In this table, m is mass of the hypothetical terrestrial exoplanet, e is 
eccentricity, i or b is the orbital inclination measured from the plane perpendicular to the sky, 
Q is the longitude of ascending node, u is the argument of pericenter, and M is the mean 
anomaly. Since our timing precisions are at the level of a signal from a terrestrial planet, we 
limited the mass coverage to the "super-Earth" range of 1 - 10 M e planets. Our semi-major 
axis coverage spans the 2:1 resonance, with one choice directly at the commensurability. All 
other parameters were varied so as to fully cover all possible architectures, but still to keep 
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the entire program tractable. Each simulation required > 1 hour of CPU time, or 
CPU hours total. 



7 million 



For the N-body simulations, we used the symplectic integrator code in the HNBody 
package ( IRauch fc Hamilton! 120021 n. which includes general relativity. We integrated each 
trial for 10 6 orbits of XO-2b (^ 72 00 years), which is eno ugh to identify > 95% of unstable 
orbits for other exoplanet system (IBarnes fc Quinnl 120041 ) . Each simulation could produce 
one of three outcomes: stable, unstable, or fail. Stable configurations lost no planets due to 
gravitational perturbations, while unstable ones did. Failed systems did no t conserve energy 
to be tter than 1 part in 10 4 , which is required for symplectic integrators (IBarnes fc Quinn 
2004h . 



In Figure El a sample slice through the data set is shown. We limited the visualization 
to cases with m ~ 1 M e planets with i or b = ^ = oo = 0, and M = 144°. Green squares 
designate stable configurations, red unstable, and yellow failed. Most simulations from 
the yellow bins are probably unstable, as close approaches between planets can violate the 
algorithm's underlying assumption that the gravitational force from the star is much larger 
than that from any other planet. We note that the average fractional change in energy 
(dE/E) for stable cases was 7 x 10 -9 . 

The trials shown in Fig. [6] show the that, for this value of M, the mean motion 
resonance at 0.0587 AU stabilizes the system. At all values of a and low e, the system 
is stable, as predicted by Hill stability theory (see iMarchal fc Bozisl 1 1982 : Idadmanl Il993l : 
Barnes fc Greenbergn2006l : iKopparapu fc Barnesll2010l ). However at the resonance, stability 
is likely for e < 0.7. As expected, we find that only certain values of M predict this "tongue" 



4 Publicly available at http://janus.astro.umd.edu/HNBody/ 



Table 9: Range of Parameter Space for Stability Calculations 



Parameter Unit Min Max 



N 



Mass 


M e 


1 


10 


1 


10 


a 


AU 


0.055 


0.062 


0.0016 


5 


e 


- 





0.8 


0.1 


9 


^orb 


o 


0.001 


70 


10 


8 


Q 


o 





324 


36 


10 


00 


o 





324 


36 


10 


M 


o 
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36 


10 
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Stability Diagram XO-2 
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Fig. 6. — Contours of TTVs for orbital configurations of a 1M companion to XO-2b, close 
to the 2:1 resonance. The color-coded squares note the stability of N-body simulations, with 
yellow (light grey) being unstable, white being failed, and green (dark grey) being stable 
simulations. The color version of this image is available in the online version. 



of stability. 

The contour lines in Fig. [6] show the values of the TTV signal in seconds. For the broad 
stable region at e < 0.1, the signal is at or below the detection limit. For larger e, especially 
in the resonance, the value can be much larger. For the bin at a = 0.0587 AU, and e = 
0.7, the signal magnitude is close to 2377 sec, or nearly 40 minutes. Clearly, the APOSTLE 
project could have detected an Earth-mass companion if it were in a favorable orbit. We do 
note, however, that an Earth-like planet with large eccentricity is likely to be rapidly tidally 



circu larized due to the close proximity (IRasio et al .11 19961 : 1 Jackson et al.ll2008l : iBarnes et al 
2010), although some eccentricity would be maintained by the resonance. 
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As noted in Table 151 (§ 14.21) the median TTV errors obtained from APOSTLE lightcurves 
of XO-2b were ~20 sec. Of the 3.6 million parameter configurations, ^1.14 million were 
found to be stable, and approximately 51% of these stable configurations produce TTVs 
on the order of 20 sec or less and thus would not be easily detected by APOSTLE or 
other surveys with similar capability. The vast amount of data generated by the stability 
simulations cannot be adequetly discussed in this text, hence we provide the results from all 
3.6 Million simulations as part of an online-only table; a segment of which is displayed in 
Table [TU1 The first column 'Sim. ID' denotes the simulation number and columns numbered 
(2)-(8) are the input configurations for a given simulation. The range of these simulation 
configurations were summarized earilier in this section and in Table [91 Columns (9)-(13) are 
the simulation statistics. Column (9) indicates the stability outcome of a given simulation, 
with or 1 denoting unstable and stable simulations respectively. Column (10) lists the 
energy conservation of a given simulaions as noted by the maximum fractional change in 
energy over a given step (dE/E). Columns (11), (12) and (13) list the standard deviations 
in the transit timing variations, transit duration variations and the variations in the impact 
parameter (6) respectively. 



6. Conclusions 



I. Photometric Precision: The XO-2 system was observed over a variety of observing 
conditions over a period of 3 years between 2008 2011. The best photometric precisions 
we obtained with our /-band and /-band observations were 405 and 354 ppm respectively 
(see Tabled]). 

II. XO-2b System Parameters: Our analysis of the 10 transit lightcurves yielded esti- 
mates of system parameters that agree with measurements presented by other studies (see 
Figure [3]). We were able to improve the constraint s on a/R+ and p + by ^40% compared to 



the previous best measurements from the ground (ISing et al.ll201ll ). The measurements are 
presented in Table [3 see the TAP values. We could not get the Markov chains to converge 
while fitting for stellar limb- darkening parameters, echoing the results from previous studies. 

III. Search for Transit Depth Variations: Our Multi-Depth fits show some variations in 
the transit depth over transit epoch (see Figure [Hand Table [6j) for the 6 /-band lightcurves. 
Since 3 out of the 6 lightcurves were not fully sampled we cannot confidently assert real 
variability in the data. The seen variations could be shortcomings of the transit model's 



Table 10: Data on the Stability Simulations of the XO-2 system 



Sim. ID 


Mass 


a 


e i Q oo 


M 


Stability 


Energy Cons. 


std(TTV) 


std(TDV) 


std(TBV) 


- 


M e 


AU 


o o o 


o 


- 


- 


seconds 


days 


- 


(1) 


(2) 


(3) 


(4) (5) (6) (7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 






1.0 


0.0555 


0.0 





0.0 








1 


2.370e-09 


2.040e+00 


2.656e-06 


2.066e-08 




1 


1.0 


0.0555 


0.0 





0.0 





36 


1 


2.630e-09 


2.062e+00 


2.633e-06 


2.000e-08 




2 


1.0 


0.0555 


0.0 





0.0 





72 


1 


2.610e-09 


2.077e+00 


2.710e-06 


2.179e-08 




3 


1.0 


0.0555 


0.0 





0.0 





108 


1 


2.050e-09 


2.068e+00 


2.639e-06 


2.165e-08 




4 


1.0 


0.0555 


0.0 





0.0 





144 


1 


1.980e-09 


2.027e+00 


2.684e-06 


2.102e-08 




5 


1.0 


0.0555 


0.0 





0.0 





180 


1 


9.120e-10 


2.029e+00 


2.662e-06 


1.983e-08 




6 


1.0 


0.0555 


0.0 





0.0 





216 


1 


2.210e-09 


2.017e+00 


2.680e-06 


2.031e-08 


to 


7 


1.0 


0.0555 


0.0 





0.0 





252 


1 


2.300e-09 


2.050e+00 


2.677e-06 


2.140e-08 


oo 


8 


1.0 


0.0555 


0.0 





0.0 





288 


1 


2.430e-09 


2.051e+00 


2.673e-06 


2.009e-08 




9 


1.0 


0.0555 


0.0 





0.0 





324 


1 


2.610e-09 


2.031e+00 


2.665e-06 


1.958e-08 





*The data are presented in their entirety as an online-only table. 

(I) Simulation ID, (2) Mass of perturber, (3) Semi-major axis, (4) Eccentricity, (5) Inclination, (6) Longitude of the ascendening node, 
(7) the Argument of Pericenter, (8) the Mean Anomaly, (9) Stability Outcome, (10) Energy Conservation dE/E, 

(II) Standard Deviation of TTV - Transit Timing Variations 

(12) Standard Deviation of TDV - Transit Duration Variations 

(13) Standard Deviation of TBV - Transit Impact Parameter (b) Variations 
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ability to fit a set of lightcurves with both complete and incomplete data. We do not see 
similar variations in the r'-band, where one may expect spot-modulated variability to appear, 
due to the greater spot-to-star contrast in the /-band. There are no known reports of stellar 
activity on XO-2, hence the /-band results are consistent with this fact. 



IV. Search for Transit Timing Variations: The XO-2b dataset contains lightcurves with 
some of the best photometric precisions achieved with ground-based observations. Since 
photometric precision directly translates to transit timing precision we are able to report 
timing precisions as low as 12 sec (after red-noise analysis, see Table [8]). We were unable to 
detect significant timing deviations for XO-2b in our data. The linear fit was not robust, 
with x 2 1 DOF being significantly greater than 1, indicating large scatter around the linear 
ephemeris fit. The transit t i mes d erived from the APOSTLE lightcurves using the red-noise 
analysis of ICarter fc Winnl (120091 ) resulted in more conservative errors than those derived 
using TMCMC. A linear ephemeris is consistent with the transit time measurements reported 
from the TAP analysis. The overall variation in the O-C values from our rednoise analysis 
was on the order of ~ 39 sec. From a sinusoidal fit to APOSTLE's O-C data, we obtained a 
TTV amplitude of ^ 27.4 sec. However, checking the goodness-of-fit of sinusoids to random 
rearrangements of the data show that the detected periodicity unlikely to be real. The 
resulting amplitude rules out planetary companions more massive than 0.2 M near the 2:1 
MMR, and larger companions near higher order resonances. 

We conclude that the set of transit times published in the literature for XO-2b and 
other transiting systems in general, are not suited for transit timing analysis. Lacking red- 
noise analysis leads to underestimated timing errors and may lead to premature reporting 
of timing variations. A proper analysis of transit times would need a simultaneous analysis 
of transit lightcurves using a transit mo del which is, (1) su ited for Bayesian inference {i.e. 
with a fairly uncorrelated parameter set, ICarter et al.ll2008f) and (2) a transit model which 
can adequately account for red-noise in the data (like TAP, iGazak et al.l 120111 ). Using large 
data sets of transit lightcurves may be inefficient due to the slowness of Markov chains with 
the addition of model parameters. In addition to the d evelopment of more detailed models 
the utilization of fast Markov chain algorithms (for e.g. iForeman-Mackey et al.ll2012l ) is also 
recommended. 

Our lack of detection of TTVs in the data for the Hot- Jupiter XO-2b is also con- 
sistent with Kepler ' s findings that (1) Hot-Jupiters tend to lack other planetary siblings 
( ILatham et al.ll201ll : ISteffen et al.ll2012l ) and (2) members of multi-plane t systems with short 



period planets (Perio d < 10 days) are more likely to be Hot-Neptunes ( ILatham et al 
Lissauer et al.ll2011bl ). 
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V. Dynamical Study: We ran 3.6 million N-body simulations of possible multi-planet 
configurations near the 2:1 resonance in order to test for (1) orbital stability and (2) de- 
tect ability of TTV signals. We varied several properties of the hypothetical companion to 
XO-2b (see Table [9]), to look for patterns in the resulting stability and transit times of the 
simulations over parameter space. Of the several stable configurations we find that ~51% 
of the simulations would display TTVs weaker than the precision limits of our survey and 
hence remain undetectable. The entire table of simulation statistics are presented online (see 
Table HUj). 
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